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I.  INTRODUCTION 


This  program  is  designed  to  simulate  the  endgame  phase 
between  a  missile  and  an  air  target  for  the  purpose  of 
illustrating  the  dynamics  of  the  encounter.  This  program, 
called  ENDGAME,  enables  the  user  to  determine  the  miss 
distance,  the  fragment  spray  zone,  and  the  best  point  for 
detonation  for  a  given  ordnance  package,  target  and 
missile  velocity  and  flight  path. 

This  program,  written  in  the  C  language,  and  using  the 
HALO  graphics  subroutine  packages,  provides  the  basic 
framework  from  which  further  simulations  of  increased 
complexity  and  sophistication  can  be  easily  implemented. 

To  establish  the  proper  perspective,  a  brief  outline  is 
given  as  follows;  the  following  chapter  gives  the 
theoretical  background  for  the  endgame.  This  is  followed 
in  Chapters  III  and  IV  by  a  construction  of  the  program  in 
the  C  language  with  comments  on  graphical  features. 
Finally  conclusions  are  drawn  in  Chapter  V,  and 
recommendations  are  made  for  future  improvements .  The 
program  variables  are  described  in  Appendix  A.  The  program 
listing  is  given  in  Appendix  B,  and  a  user's  manual  is 
given  in  Appendix  C. 


For  a  given  encounter  between  a  missile  and  a  target, 
the  guidance  system  can  use  any  one  of  several  methods  to 
navigate  the  missile  along  the  flight  path  from  launch  to 
the  intercept  with  the  target.  The  miss  distance  is 
basically  the  minimum  distance  between  the  missile  and  its 
target  during  the  encounter.  It  is  used  for  determining  the 
probability  that  a  fragment  from  the  detonation  of  the 
warhead  hits  the  target.  The  miss  distance  can  be 
evaluated  based  upon  the  encounter  conditions  shown  in 
Figure  2-1  and  assuming  the  velocity  of  both  the  missile 
and  the  target  are  constant. 


Figure  2-1.  Encounter  Conditions  in  Grobal  Coordinates 


To  evaluate  the  miss  distance,  grobal  coordinates  are 
chosen  for  convenience  and  proper  understanding.  See 
Figure  2-1. 

At  time  r 

T  +  t-  lf+m  +  s  (1.1) 

where 

T  -  the  vector  from  the  origin  to  the  target  at  r* 0 

t  *  the  vector  from  the  target  at  r=0  to  the  target 

at  T«r 

M  *  the  vector  from  the  origin  to  the  missile  at  r*=0 

m  »  the  vector  from  the  missile  at  r= 0  to  the  missile 

at  rar 

♦ 

s  *  the  vector  from  m  to  t  (the  seperation  vector) 
Since  the  target  velocity,  Vt,  is  assumed  to  be  constant 


t  -  Vt  ri  (1.2) 

And  since  the  missile  velocity,  Vffl,  is  also  assumed  to  be 
constant 

-* 

m  -  V_  (cos0  I  +  sin0  H)  ^  (1*3) 

thus, 

s  -  (T  -  M)  +  (t  -  m) 

-  [(Tx  -  Mjj)  +  Vt  r  -  Vm  COSe  r]  i  (1.4) 

+  [(Ty  -  My)  -  r  Vm  sin0]  S 


where  the  subscripts  x  and  y  denotes  the  x  and  y  components 


For  min  | s | 


d|s| 

-  =  0  (1.5) 

dr 

Solving  equation  (1.5)  for  r  using  equation  (1.4)  gives 


Sx  (Vm  cosff  “  Vt)  +  Sv  Vm  s^n& 

- ^ - 1 - - -  (1.6) 

(vm  cos#  -  vt)2  +  (Vm  sin*)2 


Define 


and 


Sx  -  Tx  -  Mx 

Sy  —  Ty  “  My 


'mtx 


V, 


mty 


Vm  cos* 
Vm  sin* 


-  V* 


Eq  (1.6)  can  be  written  as 


sx  vmtx  +  sy  vmty 


mtx 


+  V, 


mty 


(1.7) 


and  eq  (1.4)  can  be  given  in  the  form 


S  -  Csx  '  vmtx  '1  1  +  Csy  *  vmty  ')  5 


(1.8) 


Equation  (1.8)  can  be  used  to  determine  whether  the 
given  encounter  is  an  Early  Bird  or  a  Late  Bird  situation. 
In  the  Early  Bird  case,  the  missile  passes  in  front  of  the 
target,  whereas  it  passes  behind  the  target  in  the  Late 
Bird  case. 
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B.  THE  FRAGMENT  SPRAY  ZONE 


When  a  warhead  is  detonated  near  an  aircraft,  tr. 
fragments  are  usually  ejected  uniformly  and  begin 
propagate  outward  in  a  divergent  spherical  pattern  a" 
velocity  that  is  the  vector  sum  of  the  initial  fragme.- 
velocity  from  a  static  warhead  detonation  and  the  miss., 
velocity. [Ref . 1] 

If  the  fragments  are  assumed  to  have  uniform  ve.cc.t 
and  to  be  uniformly  spread  over  a  spherical  segment,  tr. 
fragment  spray  density  at  the  distance  R  is  giver,  by 

N 

p  -  - - - ,1.9 

2 jtR^  ( cos  -  COS  $2) 


where 

N  -  the  total  number  of  fragments  j.n  the  warhead 


The  equation  defining  the  angles  ^  and  <p2  is 


tan 


-1 


vm  sins  +  VQ  sin (  9  +  aL) 


L  Vm  cosd  +  VQ  cos( 
1/2 


—— i- - J  -< 

+  <*i)  ±  Vt  i 


(1.10) 


where 


Vjjj  -  the  speed  of  the  missile 
6  ■  the  elevation  angle  of  the  missile 

VQ  ■  the  average  fragment  speed  with  respect  to  a 
stationary  warhead 


q-j_  ■  the  static  leading  fragment  spray  angle 
a2  ■  the  static  trailing  fragment  spray  angle 
-  the  target  horizontal  speed 

C.  THE  BEST  POINT  FOR  DETONATION 

The  fragment  spray  density  from  the  detonation  point  is 
determined  by  the  given  detonation  distance.  At  this 
point,  the  following  problem  is  encountered.  Where  is  the 
best  point  for  detonation  ?  The  specific  answer  to  that 
question  requires  an  understanding  that  there  are  four 
different  fragment  spray  zones  created  by  the  encounter 
conditions  and  the  warhead  parameters  as  shown  in  Figure 
2-2. 


Figure  2-2.  The  Fragment  spray  zone 
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In  region  1  and  3,  part  of  the  fragment  spray  zone  hits 
part  of  the  target.  In  region  2,  the  entire  fragment  spray 
zone  hits  part  of  the  target,  and  in  region  4,  part  of  the 
fragment  spray  zone  hits  all  of  the  target. 


III. 
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III.  PROGRAMMING  "ENDGAME" 

This  chapter  contains  the  programming  description  of 

the  implementation  of  the  endgame  equations  using  the 

Lattice  C  compiler.  Even  though  several  compilers  are 

available  these  days,  the  Lattice  C  compiler  is  used  since 

it  is  a  portable  compiler  for  the  high  level  programming 

language  called  C. 

1 

A.  PREPARING  SOURCE  FILE 

I 

The  source  file  was  written  by  using  the  text  line 

I 

editor  EDLIN  in  IBM-DOS  version  3.2  to  create,  change,  and 

display  the  source  file.  All  the  external  variables  are 

declared  outside  of  the  module,  including  other  standard 

! 

i 

* 

1 

< 

i 

» 

source  files.  The  driver  program,  main() ,  is  the  first 

routine  to  be  executed.  This  is  followed  next  by  several 

subroutine  programs,  sub-subroutines,  and  so  on.  Since  the 

program  is  so  massive  that  it  is  a  little  hard  to  read  the 

whole  program  thouroughly  at  a  glance,  each  of  the 

subroutine  functions  are  entitled.  See  the  program  listing 

in  Appendix  B  for  further  particulars. 

i 

B.  COMPILING 

The  Lattice  C  compiler  version  3.00  [Ref. 3]  is  used 

S 

i 

i 

• 

for  compiling  the  source  file  into  the  object  file.  The  LC 

1 

m 

i 

1 

i 

i 

. 

command  is  used  with  several  options:  -w  option  to  shut  off 
warning  messages  generated  for  return  statements,  and  -k2 
option  to  specify  the  generate  code  for  80286.  Another 
option  chosen  is  -s  for  naming  program  segments  to  provide 
more  availability  on  RAM  working  memory.  Otherwise  the 
storage  requirements  exceed  the  available  memory. 

This  was  the  most  difficult  part  of  the  author’s  thesis 
work:  to  find  the  proper  options  that  would  work  with  the  P 
model  as  the  program  got  bigger  and  bigger.  Refer  to  the 
Lattice  C  compiler  manual  for  more  detailed  information. 

C.  LINKING 

When  linking  the  object  file  into  the  program  file,  the 
LINK  command  in  MS-DOS  version  3.0  was  used.  Since  most 
parts  of  the  program  call  the  graphics  subroutine  functions 
in  HALO  to  draw  the  graphics  on  the  screen,  a  special 
format  for  linking  the  program  with  the  HALO  library  is 
required.  With  program  model  (version  3.0)  supported  by 
the  full  available  working  memory,  the  following  command  is 
used: 

LINK  CP+ENDGAME+HALODVP3 , ENDGAME , , LCMP+LCP+HALOP3 
Further  details  on  each  step  is  discussed  in  Section  4 
of  the  HALO  manual  and  can  be  refered  to  by  looking  into 
the  LINK  command  in  the  MS-DOS  manual. 


As  mentioned  in  the  Introduction,  most  of  the  author's 
efforts  were  devoted  to  creating  screens  that  illustrate 
the  endgame  geometry  using  colorful  graphics  to  obtain 
desirable  results  from  the  given  warhead  and  encounter 
data  input.  In  ENDGAME,  the  useful  information  is 
presented  on  each  screen  with  high  resolution,  is  processed 
and  displayed  graphically  on  the  screen,  simulating  each 
endgame  phase  and  showing  what  is  happening. 

All  the  screens  were  created  using  HALO  in  the  graphics 
mode,  except  the  option  screen  which  shows  the  sequential 
questions  and  input  data.  The  HALO  library  used  in  this 
work,  version  2.26,  is  a  collection  of  high  performance 
subroutines  which  allow  the  application  programmer  '  to 
develope  sophisticated  computer  graphics  programs. 

ENDGAME  uses  world  coordinates  throughout  the  program. 
The  world  coordinates  are  mapped  to  device  coordinates  on  a 
specified  range  from  0  to  640  on  x-axis  and  from  o  to  350 
on  y-axis  so  as  to  locate  each  point  within  the  resolution 
of  the  EGA (Enhanced  Graphics  Adaptor)  display  device.  The 
upper  left  corner  of  the  screen  is  at  coordinates  (0,0)  and 
the  lower  right  corner  is  at  coordinates  (639,349)  for 
convenience.  To  draw  objects,  the  coordinates  and 
position  of  the  objects  have  to  be  scaled  by  the  aspect 


i»JA 


ratio.  Aspect  is  expressed  as  a  ratio  of  width  to  height. 
The  value  1.439  was  emperically  chosen  to  make  the  range 
increment  on  the  horizontal  axis  be  the  same  as  that  on 
the  vertical  axis. 

Another  factor  considered  was  trans_f actor  [Ref.2],  As 
mentioned  in  the  user's  manual  in  Appendix  C,  some 
particular  situations  had  to  be  excluded  in  the  endgame. 
The  user  should  exclude  the  encounter  conditions  where  the 
difference  in  relative  position  between  the  target  and  the 
missile  is  larger  than  1,000  ft.  Hence,  the  difference 
always  lies  within  1,000  ft  and  the  graphic  resolution  is 
taken  account  into  by  trans_f  actor.  The  value  of 
trans_f actor  such  as  0.2,  0.3,  etc.,  is  emperically  chosen 
based  upon  the  value  that  gives  good  graphic  resolution 
on  the  screen. 


A.  CREATED  SCREENS 

All  the  screens  created  by  HALO  are  in  the  graphics 
mode.  There  are  three  different  screens  which  allow 
students  to  comprehend  the  dynamic  situation  between  the 
target  and  the  missile. 


When  ENDGAME  is  executed,  the  opening  screen  is 
displayed,  welcoming  the  student  to  the  endgame  assessment 
as  shown  in  Figure  4-1.  Hitting  any  key  allows  the  student 
to  proceed  to  the  first  step. 


cv:\* 
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After  inputting  the  data  requested  by  the  option 
screen,  the  student  will  see  the  miss  distance  screen 
(Figure  4-2) .  This  screen  tells  the  student  what  the  miss 
distance  is  numerically  and  shows  how  it  looks 
vectorically  with  the  target  and  the  missile  vectors  shown 
on  the  screen.  It  also  prints  and  shows  whether  the 
encounter  situation  is  an  Early  Bird  or  a  Late  Bird.  Refer 
to  the  Appendix  A  for  the  other  variables  shown  on  the 
screen. 

3 .  The  Fragment  Sorav  Zone  Screen 

This  screen  (Figure  4-3)  shows  the  fragment  spray 
zone  with  respect  to  the  target.  Unfortunately,  the  screen 
doesn't  show  the  different  fragment  vectors  due  to  leading 
and  trailing  fragment  spray  angle  with  different  colors. 
However,  the  reader  should  be  able  to  visuallize  the 
fragment  pattern.  Also  refer  to  the  Appendix  A  for  other 
variables. 

4.  The  6?st  pgla.t...f<ar  Detection.  Screen 

The  student  will  see  where  the  best  point  for 
detonation  is  on  this  screen  (Figure  4-4) .  It  shows  the 
zones  which  are  appropriate  for  the  detonation  of  the 
warhead  for  the  given  encounter  situation  considering  the 
combination  of  the  dynamic  fragment  speed  with  respect  to 
the  target  and  the  missile  velocity  vector.  From  this 
screen,  the  student  can  go  back  to  see  the  previous  screens 


•imply  by  following  the  instructions  on  the  bottom  of  the 
screen.  Finally,  the  end  screen  will  appear  if  the  student 
follows  the  appropriate  instruction  on  the  bottom  of  the 
screen. 

B.  OPTION  SCREENS 

Several  questions  appear  on  the  screen  prompting  the 
student  to  enter  the  necessary  input  data.  This  data  is 
used  to  calculate  the  miss  distance  and  the  best  point  for 
detonation.  These  screens  are  not  in  the  graphics  mode,  but 
are  in  standard  text  mode.  Warning  or  error  messages  might 
appear  if  the  student  inputs  unreasonable  data.  The 
limitations  for  this  endgame  program  are  described  in  the 
user's  manual  which  is  given  in  Appendix  C. 
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Figure  4-3.  The  Fragment  Spray  Zone  Screen 


This  computer  graphics  simulation  program  for  the 
endgame  between  a  missile  and  an  air  target  was  written  in 
the  C  language  using  the  Lattice  C  compiler  running  on  MS- 
DOS.  The  program  is  divided  into  three  major  categories: 
miss  distance,  fragment  spray  zone,  and  best  point  for 
detonation.  It  was  intended  to  visually  demonstrate  the 
dynamic  evaluation  for  the  terminal  events  in  the  encounter 
between  an  aircraft  and  a  high  explosive  proximity  fuzed 
warhead  on  the  missile. 

It  is  hoped,  finally,  that  students  who  are  more 
interested  in  this  field  will  continue  to  develop  the 
present  ENDGAME  program  with  an  anxiety  and  an  eagerness. 
Especially,  a  more  realistic  endeavor  to  expand  the  variety 
of  options  within  the  program  to  account  for  any  encounter 
situation  is  strongly  recommended.  For  example,  adding  a 
target  acceleration  to  the  target  motion  would  allow  more 
realism.  From  a  programming  point  of  view,  a  solution  to 
the  memory  expansion  problem  is  also  desirable.  Hopefully, 
the  probability  of  kill  will  be  developed. in  future  study. 


APPENDIX  A 


PROGRAM  VARIABLES 

The  variables  used  in  endgame  computer  program  are 


listed 

as  follows; 

alphal 

alpha2 

[degree] 

[degree] 

static  leading  fragment  spray  angle 
with  respect  to  missile  axis 
static  trailing  fragment  spray  angle 
with  respect  to  missile  axis 

gamma 1 

gamma 2 

[degree] 

[degree] 

angle  from  horizontal  to  dynamic 
leading  fragment  spray 
angle  from  horizontal  to  dynamic 
trailing  fragment  spray 

ht 

[degree] 

target  heading 

length 

[ft] 

target  presented  length 

minsx  [ft] 
minsy  [ft] 
miss_distance  [ft] 
mx  [ft] 

my  [ft] 

x  component  of  miss  distance  vector 
y  component  of  miss  distance  vector 
miss  distance  or  closest  point 
x  component  of  missile  position 
y  component  of  missile  position 

num 

total  number  of  fragments  in  warhead 

phione 

phi two 

[degree] 

[degree] 

dynamic  leading  fragment  spray  angle 
with  respect  to  missile  axis 
dynamic  trailing  fragment  spray  angle 
with  respect  to  missile  axis 

r_det 

rho 

[ft] 

[frag/ ft2] 

detonation  distance 
fragment  spray  density 

sx 

sy 

[ft] 

[ft] 

x  component  of  difference  between 
missile  and  target  position 
y  component  of  difference  between 

if  'll  ’** 


missile  and  target  position 


tau 

[sec] 

minimum  time  for  miss  distance 

theata 

[degree] 

elevation  angle  of  missile 

tx 

[ft] 

x  component  of  target  position 

ty 

[ft] 

y  component  of  target  position 

vftl 

[ft/sec] 

dynamic  leading  fragment  speed  with 
respect  to  target 

Vft2 

[ft/sec] 

dynamic  trailing  fragment  speed  with 
respect  to  target 

vftxl 

[ft/sec] 

x  component  of  vftl 

VftX2 

[ft/sec] 

x  component  of  vft2 

vftyl 

[ft/sec] 

y  component  of  vftl 

vfty2 

[ft/sec] 

y  component  of  vft2 

vm 

[ft/ sec] 

speed  of  missile 

vmtx 

[ft/ sec] 

x  component  of  missile  speed  with 
respect  to  target 

vmty 

[ft/sec] 

y  component  of  missile  speed  with 
respect  to  target 

VO 

[ft/ sec] 

average  fragment  speed  with  respect 
to  stationary  warhead 

vt 

[ft/sec] 

target  horizontal  speed 
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/*****************************************************/ 
/******  define  the  all  external  variables  *******/ 

/*****************************************************/ 


/******  include  another  source  file  ******/ 

# include  "stdio.h" 

# include  "stdlib.h" 

♦include  "math.h" 

/**********  image  file  memorized  *********/ 

static  char  filel[ ]={ "miss. pic" } ; 
static  char  file2 [ ]={ "prob.pic" } ; 
static  char  file3 []«{ "best. pic" } ; 


/**********  position  variables  ***********/ 

float  xl=0.0,yl=*0.0,x2=640.0,y2=350.0;  /*world  coordinate*/ 
float  gcx,gcy  ;  /*  graphic  cursor  position  */ 

float  gcxl , gcx2 , gcyl , gcy2 ;  /*  for  varing  graphic  cursor  */ 

float  tcx,tcy;  /*  text  cursor  position  */ 

float  dx,dy;  /*  cursor  increment  */ 

float  bxl , byl , bx2 , by2  ;  /*  box  creation  */ 

/***  variables  for  adjusting  to  screen  ***/ 


float  trans_factor  ; 

float  dist_factor,dist_l,dist_2 ? 

/************  color  definition  ***********/ 


int 

black  =0; 

int 

blue  *  1  ; 

int 

green  =  2 

• 

t 

int 

Cyan  ■  3  ? 

int 

white  *  7 

• 

9 

int 

liblue  -  9 

• 

9 

int 

ligreen  - 

10 

int 

licyan  - 

11 

int 

lired  - 

12 

int 

lipurple  - 

13 

int 

liyellow  - 

14 

/***********  system  variables  ************/ 

float  aspect=0.7249;  /*  ratio  of  width  to  height  */ 
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int  mode  -  4;  /*  enhanced  graphics  adapter  mode  */ 

int  page-2;  /*  available  graphics  page  for  system  */ 


/************  drawing  variables  *******■***/ 


float  angl,ang2; 
int  width; 
int  style; 
float  radius; 


/*  starting  and  ending  angle  of  arc  */ 

/*  line  width  */ 
/*  line  style  */ 
/*  radius  of  arc  in  world  coordinates  */ 


char  s[15] ; 
int  c; 

int  h,w,p,m; 

float  itx,ity; 
float  imx , imy ; 
int  count; 


/*********  variables  for  program  *********/ 

float  tx,ty  ; 

float  mx,my  ; 

float  sx,sy  ; 

float  vt,vm  ; 

float  vo  ; 

float  ht  ; 

float  theata; 

float  tau; 

float  vmtx,vmty  ; 

float  minsx,minsy  ; 

float  m_dist; 

float  vftxl , vftyl ; 

float  vftx2,vfty2; 

float  gamma l,gamma2; 

float  alphal,alpha2 ; 

float  vftl,vft2; 

float  phione,phitwo  ; 

float  r_det,num, length, rho; 
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z*********************************************************y 

/**  main  function  **/ 

/*********************************************************y 


SP 


W  W1  ^  Itfiim  r^;wv«TA''a«-mTTX"j»-\.»ron-3tr 


/★ft*******************************************************/ 

/**  SUBROUTINE:  OPENING ()  **/ 

/********************************************************* / 

openinq ( ) 

{ 

varl3 ( ) ; 
varl4 ( ) ? 

setcolor (&lired) ; 
gcx*3  0.0, gcy«3  20.0; 
movabs ( &gcx, &gcy ) ; 
gcx«320.0; 
lnabs ( &gcx, &gcy ) ; 

settextclr (&lired,&black) ; 
tcx*340 . 0 , tcy*320 . 0 ; 
movtcurabs ( &tcx, &tcy) ; 

text ("Subject  :  AIR  DEFENCE  LETHALITY"); 

setcolor (Slipurple) ; 
gcx*4  0.0, gcy»3 10.0; 
movabs ( &gcx, &gcy) ; 
gcx«330.0; 
lnabs ( &gcx, &gcy ) ; 

settextclr (Slipurple, Sblack) ; 
tcx-350 . 0 , tcy=3 10 . 0 ; 
movtcurabs (&tcx, &tcy) ; 
text ("Professor  :  R.  E.  Ball"); 

h«4 , w»4 ; 

settext (&h, &w, &p, &m) ; 
settextclr ( figreen, iblack) ; 
tcx*2 10 . 0 , tcy-2 2 0 . o ; 
movtcurabs ( &tcx, &tcy) ; 
text ( "WELCOME" ) ; 

h-3,w»3; 

settext ( &h ,  &w ,  &p ,  &m ) ; 
settextclr (fcligreen, Sblack) ; 
tcx-299 . 0 , tcy-185 . 0 ; 
movtcurabs ( &tcx, &tcy) ; 
text ("to"); 

h-6,w«6; 

settext (&h,&w,&p,&m) ; 
settextclr (filicyan, Sblack) ; 
tcx«155 . 0 , tcy-107 . 0 ; 
movtcurabs (&tcx,&tcy) ; 
text ( "ENDGAME" ) ; 
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h»2 , w-2 ; 

settext (&h, &w, ip, &m) ; 

settextclr(&liyellow, Sblack) ; 

tcx-123 . 0 , tcy-55 . 0 ; 

movtcurabs ( &tcx, &tcy) ; 

text ( "NAVAL  POSTGRADUATE  SCHOOL"); 

bottom ( ) ; 

varl2 ( ) ; 

) 


/*********************************************************/ 

/**  SUBROUTINE:  0PTI0N1()  **/ 

/*********************************************************/ 

optionl ( ) 

{ 

labl: 

printf("\n  Initial  x-position  of  target  [ft]  ? 
scanf ("%f" ,&tx) ; 
if((tx<0.0) | | (tx>100000.0) ) 

{ 

printf("\n  Unreal!  Assume  that  it  be  above  ground  and 
real . " )  ; 

printf("\n  Type  again."); 
goto  labl; 

} 

else 

{ 

lab2 : 

printf("\n  Initial  y-position  of  target  [ft]  ?  :"); 

scanf ("%f",&ty) ; 

if ( (ty<0.0) | | (ty>100000.0) ) 

{ 

printf ("\n  Unreal!  Assume  that  it  be  above  ground  and 
real . " )  ; 

printf ("\n  Type  again."); 
goto  lab2; 

} 

else 

{ 

lab3 : 

printf ("\n  Target  speed  [ft/sec]  ?  :"); 

scanf ("%f",&vt) ; 

if ( (vt<0.0) | | (vt>2000) ) 

{ 

printf ("\n  Unreal!  Type  again."); 
goto  lab3 ; 

} 

else 
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{ 

lab4: 

printf("\n  Target  heading  [degree]  ?  (  type  0  or  180  ) 

.  N  )  . 

scanf  ("%f",&ht)  ; 
if ((htl-0)&&(hti-180) ) 

{ 

printf(H\n  Assume  level  flightl  Type  again."); 
goto  lab4; 

} 

else 

{ 

lab5: 

printf("\n  Initial  x-position  of  missile  [ft]  ? 

scanf ("%f",&mx) ; 

if ( (mx<0. 0) | | (mx>100000. 0) ) 

{ 

printf(M\n  Unreal!  Assume  that  it  be  above  ground  and 
real  . " )  ; 

printf("\n  Type  again."); 
goto  lab5; 

} 

else 

{ 

lab6: 

printf("\n  Initial  y-position  of  missile  [ft]  ?  :"); 

scanf ("%f" ,&my) ; 

if ( (my<0.0) | | (my>100000 . 0) ) 

{ 

printf("\n  Unreal!  Assume  that  it  be  above  ground  and 
real . " ) ; 

printf("\n  Type  again."); 
goto  lab6; 

} 

else  if (( (abs(tx-mx) <100.0) && (abs (ty-my) <100.0) )  || 

( (abs(tx-mx) >1000.0) && (abs (ty-my) >1000.0) ) ) 

{ 

printf("\n  Unproper  input  data  for  end  game!  Start 
again. ") ; 

goto  labl; 

) 

else 

{ 

lab  7 : 

printf("\n  Missile  speed  [ft/sec]  ?  :"); 

scanf ("%f" , &vm) ; 

if ( (vm<vt) | | (vm>6000. 0) ) 

{ 

printf("\n  Unreal!  Missile  speed  must  be  greater  than 
target  speed . " ) ; 

printf("\n  Make  it  real!  Type  again."); 


goto  lab7; 

} 

•lse 

{ 

lab8: 

printf("\n  Elevation  angle  of  missile  (theata)  [degree]  ? 

•  M  )  • 

scanf theata) ; 

if  ( 

( ( (ht— 0)  &&  ( ( ( (txcmx)  &&  (ty>my) )  &&  ( (180<»theata)  &&  (360>=thea 


( ( (tx<mx)  &&  (tycmy) )  &&  ( (0<*»theata)  && (theata<»180) ) ) ) ) 

I  I 

( (ht—0)  &&  ( ( ( (tx>mx)  &&  (ty>my) )  &&  ( (theata>*90)  &&  (theata<=360 

))) 

|  |  ( ( (tx>mx)  &&  (ty<my) )  &&  ( (theata>**0)  &&  (theata<*270) ) ) ) ) ) 

I  I 

( ( (ht— 180)  &&  ( ( ( (tx<mx)  &&  (ty>my)  &&  ( (theata>=*180)  &&  (theata<= 
90))) 

I  I 

( ( (tx<mx) && (ty<my) ) && ( (theata>-270) && (theata<-180) ) ) ) ) 

I  I 

( (ht*«180) && ( ( ( (tx>mx) && (ty>my) ) && ( (theata>“270) && (theata<= 
360))) 

I  I 

( ( (tx>mx) && (ty<my) ) && ( (theata>-0) && (theata<“180) ) ) ) ) ) 

)) 

{ 

printf("\n  Target  might  not  be  hit.  Put  appropriate  angle 
!iM); 

goto  lab8; 

} 

else 

l 

ht  -  ht  *  PI  /  180.0; 
theata  ■  theata  *  PI  /  180.0? 
sx  -  tx  -  mx  ; 
sy  *  ty  -  my  ? 
if  (ht— 0) 

vmtx  ■  vm  *  cos (theata)  -  vt  ; 
else 

vmtx  ■  vm  *  cos (theata)  +  vt  ; 
vmty  -  vm  *  sin (theata)  ; 
tau  »  (  sx  *  vmtx  +  sy  *  vmty  )/( 
pow(vmtx, 2 . 0) +pow(vmty, 2 . 0) ) ? 
minsx  -  sx  -  tau  *  vmtx  ; 
minsy  -  sy  -  tau  *  vmty  ? 

m_dist  -  sqrt (pow(minsx, 2 . 0)  +  pow(minsy, 2 . 0) ) ; 

ht  -  ht  *  180.0  /  PI  ; 

theata  ■  theata  *  180.0  /  PI  ? 


change: 

printf("\n  Would  you  like  to  change  these  values  ?  (  Y/N  ) 
:"); 

c-getch ( ) ; 
if  (  c— 'y') 

{ 

printf("\n  Yes.  Hit  any  key  to  change." ); 
varl2 ( ) ; 
goto  labl; 

} 

else  if(  c== ' n ' ) 

{ 

printf("\n  Hit  any  key  to  continue."); 

} 

else 

{ 

printf("\n  Type  again  correctly."); 
goto  change; 

} 

varl2 ( ) ; 

} 


//*********************************************************y' 

/**  SUBROUTINE:  DRAW1()  **/ 

/*********************************************************^ 

drawl ( ) 

{ 

top() ; 

text ("MISS  DISTANCE  [ft]  -") ; 

input 1 ( ) ; 

inputll() ; 

coord_l ( ) ; 

loc_l ( ) ; 

vectorl ( )  ; 

indexl (); 

outputl ( ) ; 

outputll ( )  ; 

outputl2 ( )  ; 

} 


/*********************************************************/ 
/**  SUBROUTINE:  OPTION2()  **/ 

/********************************************************* / 

option2 ( ) 

{ 

labll: 

printf("\n  Fragment  velocity  (Vo)  [ft/sec]  ?  :”); 

scanf (”%f” ,&vo) ; 

if  ( (vo<0.0) | | (vo>200000.0) ) 

{ 

printf (”\n  Unreal 1  Assume  that  it  be  real.”); 
printf (”\n  Type  again.”); 
goto  labll; 

} 

else 

{ 

labl2 : 

printf (”\n  Spray  angle  (alphal)  [degree]  ?  :”); 

scanf ( "%f " , Salphal) ; 

if  ( (alphaKO.O)  |  |  (alphal>l80 . 0) ) 

{ 

printf (”\n  Unreal!  Type  again.”); 
goto  labl2; 

} 

else 

{ 

labl3 : 

printf (”\n  Spray  angle  (alpha2)  [degree]  ?  :"); 

scanf (”%f",&alpha2) ; 

if  ( (alpha2<0 . 0) | | (alpha2>180 . 0) ) 

{ 

printf (”\n  Unreal!  Type  again.”); 
goto  labl3; 

} 

else 

{ 

labl4: 

printf (”\n  Number  of  fragments  (N)  ?  :"); 
scanf (”%f”,&num) ; 
if (num<0. 0) 

{ 

printf (”\n  Unreal!  Type  again.”); 
goto  labl4; 

} 

else  • 

{ 

labl5: 

printf ("\n  Detonation  distance  (R)  [ft]  ?  :”); 
scanf (”%f” , &r_det) ; 
if  (r_det<0.0) 
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{ 

printf(M\n  Unreal 1  Type  again."); 
goto  labl5; 

} 

else 

{ 

labl6: 

printf("\n  Target  presented  length  (L)  [ft]  ? 

scanf ("%f",&length) ; 

if ( (length<0.0) | | (llength>200. 0) ) 

{ 

printf("\n  Unreal 1  Type  again."); 
goto  labl6; 

} 

else 

{ 

alphal  *  alphal  *  PI  /  180.0  ; 
alpha2  *  alpha2  *  PI  /  180.0  ; 

varl5 ( ) ; 

if (ht!«o. 0) 
vt=-vt ; 
else 
vt=vt ; 

vftxl  =  vm  *  cos(theata)  +  vo  *  cos(theata  +  alphal) 

vt  ; 

vftx2  »  vm  *  cos(theata)  +  vo  *  cos(theata  +  alpha2) 

vt  ; 

vftyl  =  vm  *  sin(theata)  +  vo  *  sin(theata  +  alphal) 

vfty2  »  vm  *  sin (theata)  +  vo  *  sin(theata  +  alpha2) 

gammal  «  atan(vftyl/vftxl)  ; 

gamma2  =  atan(vfty2/vftx2)  ; 

vftl  *  sqrt (pow (vftxl, 2 . 0)  +  pow( vftyl, 2 . 0) ) ; 

vft2  »  sqrt (pow(vftx2 , 2 . 0)  +  pow(vfty2 , 2 . 0) ) ; 

varl6 ( ) ; 

phione  *  gammal  -  theata  ; 
phitwo  *  gamma2  -  theata  ; 

rhoa(num/ (2*PI* (pow(r_det , 2 . 0) ) * (cos (phione) - 
cos (phitwo) ) ) ) ; 

) 

} 

} 

} 

) 

) 
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Iabl7  s 

printf(M\n  Would  you  like  to  change  these  values  ? 
(Y/N) : M) ; 
c-getch ( ) ; 
if  (c—'y') 

{ 

printf(M\n  Yes.  Hit  any  key  to  continue. *’)  ; 
varl2 ( ) ; 
goto  labll; 

} 

else  if  (c^-'n1)? 

{ 

printf(M\n  Hit  any  key  to  continue. ") ; 

} 

else 

{ 

printf("\n  Type  again  correctly . M) ; 
goto  labll  •, 

} 

varl2 ( ) ? 

} 


/*********************************************************/ 
/**  .SUBROUTINE:  DRAW2()  **/ 

/*********************************************************/ 

draw2 ( ) 

{ 

top()  ; 

t ext ( " FRAGMENT  S  PRAY  Z ONE " ) ; 

input2 ( ) ; 

input21 () y 

loc_2 ( ) ; 

vector2 ( ) ; 

vector21() ; 

index2 ( ) ; 

output 2 ( ) ; 

output2 1 ( ) ; 

) 


/****************************************************★****/ 

/**  SUBROUTINE:  DRAW3()  **/ 

/♦♦♦♦♦♦★♦★♦♦ilr*****************************!***************#/ 

draw3 ( ) 

{ 

top() ; 


text ("BEST  POINT  for  DETONATION"); 

loc_3 ( ) ; 

vector3 ( ) ; 

index3 ( ) ; 

output 3 ( ) ; 

output31()  ; 

coord_3 ( ) ; 

repeat () ; 

} 

/*********************************************************/ 
/*  SUB-SUBROUTINE  INPUT1()  for  drawl ()  */ 

/A********************************************************/ 

input 1 ( ) 

{ 

varl4 ( ) ; 

setcolor (Swhite) ; 

bxl-476. 0 ; by 1=165 . 0 ;bx2=638 . 0 ; by 2=295 . 0 ; 
gcx»577.0;gcy=294 .0; 
boxf  ( )  ; 
varl3 ( ) ; 

settextclr (Sligreen, Sblue) ; 
tcx*540 . 0 ; tcy=282 . 0 ; 
movtcurabs ( &tcx , &tcy ) ; 
text ("DATA") ; 

settextclr (Sliyellow, &blue) ; 
tcx-510.0;tcy=269.0; 
movtcurabs ( &tcx, &tcy) ; 
text ("Target") ; 
settextclr (Swhite, Sblue) ; 
tcx=490 . 0 ; tcy=258 . 0 ; 
movtcurabs ( &tcx, &tcy) ; 
text("l.Tx  =" ) ; 
tcy-247 . 0 ; 

movtcurabs (&tcx, &tcy) ; 
text ( " 2 . Ty  ="); 
tcy=2  36.0; 

movtcurabs ( &tcx , &tcy ) ; 
text ( " 3 . Vt  ="); 
tcy»225.0; 

movtcurabs (&tcx, & tcy) ; 
text("4.Ht  -"); 
settextclr (filiyellow, Sblue) ; 
tcx-510. 0;tcy-214 .0; 
movtcurabs ( &tcx , & tcy ) ; 
text ("Missile") ; 
settextclr (Swhite, Sblue) ; 
tcx-490. 0 ;tcy=203 . 0 ; 
movtcurabs (&tcx, & tcy) ; 
text ( " 5 . Mx  -"); 
tcy- 19 2 . 0 ; 
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aovtcurabs ( &tcx , &tcy) ; 

taxt("6.Hy 

tcy*181.0? 

aovtcurabs  ( &tcx,  &tcy) ; 
text("7.Vm  «"); 
tcy-170.0; 

aovtcurabs ( &tcx , &tcy ) ; 
text ( "  8 . theata- H ) ; 

} 

/*********************************************************/ 
/*  SUB-SUBROUTINE  INPUT11()  for  drawl ()  */ 

/*********************************************************  / 


lnputll() 

{ 

settextclr(&whlte/&blue) ; 
sprlntf (s, "%-8.1fM,tx) y 
tcx-555 . 0 ; tcy-2  58.0; 
varlO ( ) ; 

sprintf (s, "%-8.1f",ty) ; 

tcy«247.0; 

varlO ( ) ? 

sprintf  (s,  "%-6. If', vt) ; 

tcy-236.0; 

varlO () ; 

sprintf (s,"%-5. If", ht) ; 

tcy»225.0; 

varlO  ()  ; 

sprintf <s,"%-8. If ",mx) ; 

tcy-203.0; 

varlO ( ) ; 

sprintf ( s , " %-8 . If " , my ) ; 

tcy»192.0; 

varlO () ; 

sprintf (s , "%-6. if" ,vm) ; 

tcy»181.0; 

varlO ( ) ? 

sprintf (s,"%-5. If", theata) ; 
tcx-580 . 0 ; tcy-170 . o ; 
varlO ( ) ? 
ht«ht*PI/l80 . 0 ; 
theata-theata*PI/l80 . 0 ; 

} 
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/*********************************************************/ 
/*  SUB-SUBROUTINE  COORD_l()  for  drawl ()  */ 

/*********************************************************/ 


coord_l ( ) 

{ 

gcx-60 . 0 ; gey- 60 . 0 ; 
aovabs ( &gcx ,  Sgey) ; 
gcx-423.5; 
lnabs (Sgex, Sgey) ; 
gcx-80.0;gcy-50.0; 
aovabs ( Sgex, Sgey ) ; 
gcy-300.0? 
lnabs ( Sgex, Sgey) ; 

for  ( count-0 ; count< 6 ; count++ ) 
{ 

gcx-gcx+40*l .439 ; gey- 6 2 . 5 ; 
aovabs ( Sgex , Sgey) ; 
gcy-60.01? 
lnabs ( Sgex , Sgey) ; 

) 

gcy-60.0? 

for  ( count»0  ?  count<6 ; count++ ) 
{ 

gcx-80 . 01 ;gcy-gcy+40 . 0 ; 

aovabs ( Sgex , Sgey) ; 

gcx-84.0; 

lnabs ( Sgex, Sgey) ; 

) 

settextclr(Sliyellow) ; 
tcx-230 . 0 ; tcy-41 . 0 ; 
aovtcurabs ( Stex, Stcy ) ; 
text ( "x  [ft]")? 
tcx<*30 . 0  ;tcy-200. 0  ? 
aovtcurabs ( Stex, Stcy) ; 
text ( "y  [ft]"); 
sattaxtclr(Sllred) ; 
tcx-3  50 . 0  ?  tcy-37 . 0 ; 
aovtcurabs ( Stex , Stcy) ; 
tsxt(" [scale-200]") ; 

) 


/A******************************************************** / 

/*  SUB-SUBROUTINE  LOC_l()  for  drawl ()  */ 

/*********************************************************/ 

loc_l() 

{ 

var3 ( ) ; 
acftl() ; 
msl ( )  ; 

} 


/ft********************************************************/ 

/*  SUB-SUBROUTINE  VECTORIQ  for  draw()  */ 

/*********************************************************y 

vectorl ( ) 

{ 

dist_factor»abs(sx) /(cos(atan(abs(sy/sx) ) ) ) ; 
dist_l«dist_factor*trans_f actor/2 . 5 ; 
dist_2«dist_factor*trans_f actor; 

/ * ** ** * ** ** *  vector  for  Vin  **********/ 

setcolor (Slicyan) ; 

dx«di3t_l*cos (theata) ;dy=dist_l*sin(theata) ; 
dx-dx*l. 439 ; 
lnrel(&dx,&dy) ; 

/***********  vector  for  Vt  **********/ 
setcolor (&whlte) ; 
if  (  ht  — i  0.0  ) 

dx“-dist_l*vt/vm ; 
else 

dx»d i  s  t_l * vt/ vm ; 
dy«0 . 0 ; 
dx*dx*1.439; 
lnrel ( &dx, &dy) ; 

/***********  vector  for  Vint  ***#******/ 

inqgcur ( &gcx, &gcy , Swhite) ; 

gcxl-gcx ; gey 1-gcy ; 

setcolor (Sliyellow) ; 

gcx»imx+8  0.0; gcy« imy+  60.0; 

movabs ( &gcx , &gcy) ; 

lnabs ( &gcxl , &gcyl) ; 

/ *********  extension  of  Vint  **********/ 
if  (  ht  —  0.0  ) 

dx»dist_2*cos (theata) -dist_2*vt/vm; 

else 
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dx-dist_2*cos (theata) +dist_2*vt/vm; 
dy-dist_2*sin (theata) ; 
dx-dx*1.439; 
style-3 ; 

setlnstyle(&style) ; 
lnrel(&dx,&dy) ; 
style-1 ; 

setlnstyle(&style) ; 

} 


/*********************************************************/ 

/*  SUB-SUBROUTINE  INDEX1()  for  drawl ()  */ 

/*********************************************************/ 

indexl ( ) 

{ 

dx-22.0;dy-0.0; 
gcx»12  0 . 0 ; gcy-3  00 . 0 ; 
tcx-150 . 0 ; tcy=296 . 0 ; 
in_vm ( ) ; 

gey— 290. 0; 
tcy-286. 0; 
in_vt ( ) ; 

gcy-280. 0; 
tcy-276.0; 
in_vmt ( ) ; 

gcy-270. 0; 
tcy-266. 0; 
in_miss ( ) ; 

} 


/*  SUB-SUBROUTINE  OUTPUTl()  for  drawl ()  */ 

/*********************************************************/ 

outputl ( ) 

•{ 

setcolor (fiwhite) ; 

bxl-476.0;byl-70. 0;bx2=638 . 0yby2-163 .0; 
gcx-577 . 0 ;gcy-162 . 0 ; 
boxf ( ) ? 


settextclr (&ligreen, &blue) ; 


tcx-535 . 0 ; tcy-152 . 0 y 
movtcurabs ( &tcx, &tcy) ; 
text ( "OUTPUT") ; 
settextclr (Swhite, &blue) ; 
tcx-490. 0ytcy»141. Oy 
movtcurabs ( &tcx, &tcy) ; 
text ("9.  Sx  -»); 
tcy-130.0; 

movtcurabs ( &tcx, &tcy) y 
text("10.Sy  «")  y 
tcy-119 . 0 ; 

movtcurabs ( &tcx, &tcy ) ; 
text ( " 11 . Vmtx  -") ; 
tcy»108 . 0 ; 

movtcurabs ( &tcx, &tcy) ; 
text ("12 .Vmty  =") ; 
tcy*97.0; 

movtcurabs ( &tcx , &tcy ) ; 
text("13.tau  «") ; 
tcy=86.0; 

movtcurabs ( &tcx, &tcy) ; 
text("14.minsx=") ; 
tcy-75.0; 

movtcurabs ( &tcxf  &tcy) ; 
text ("15.minsy=") ; 

} 


/*********************************************************/ 
/*  SUB-SUBROUTINE  0UTPUT11()  for  drawl ()  */ 

/ ****************** *********************  ******* ***********/ 

output 11 () 

{ 

sprintf (s, "%-9. If" ,sx) ; 
tcx*572 . 0 ; tcy»14 1.0; 
varlO ( ) y 

sprintf (s, "%-9. If", sy) ; 

tcy-130.0y 

varl0() y 

sprintf (s,"%-7. 2 f", vmtx) ; 

tcy»119.0; 

varlO ( ) y 

sprintf (s,"%-7. 2 f", vmty) ; 
tcy-108 . 0 ; 
varlO () y 

sprintf (s,"%-6.4f",tau) ; 

tcy-97.0; 

varlO () y 

sprintf (s,"%-7. 2 f",minsx) y 
tcy-86. Oy 
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varlO ( ) ; 

sprintf (s, "%-7.2f",minsy) ; 

tcy-75.0; 

varlO() ; 

sprintf (a, "%-7.3f”,m_dist) ; 
h-2;w-l; 

set text (&h, &w, &p, &m) ; 
tcx»360. 0;tcy«320. 0; 
varlO ( ) ; 

} 


/*********************************************************/ 
/*  SUB-SUBROUTINE  OUTPUT12()  for  drawl ()  */ 

/A********************************************************/ 

outputl2 ( ) 

{ 

/******  miss  distance  **********/ 

setcolor (Sliblue) ; 

gcx-itx+80 . 0 ;gcy«ity+60 . 0  ? 

movabs ( &gcx , &gcy ) ; 

dx«-minsx* 1.439  * trans_f actor ; 

dy*-minsy *trans_f actor ; 

lnrel (&dx, &dy) ; 

/**  determine  for  early  or  late  bird  ***/ 

settextclr(&green, &black) ; 

tcx-3  40.0; tcy«2  80.0; 

movtcurabs ( &tcx, &tcy) ; 

deltcur() ; 

if ( ( (ty>my)&&( ( (ht=*0. 0) && (atan(vmty/vmtx) <atan(sy/sx) ) ) 

| | ( (ht>3.12) &&(atan(vmty/vmtx)>atan(sy/sx) ) ) ) ) 

I  1 

( (ty<my)  && ( ( (ht— 0. 0)  &&  (atan(vmty/vmtx)  >atan(sy/sx) ) ) 

| | ( (ht>3 . 12) && (atan(vmty/vmtx) <atan(sy/sx) ) ) ) ) ) 
text ("Early  Bird  !"); 
else 

text ("Late  Bird  !"); 
gwrite(f ilel) ; 
bottom ( ) ; 
varl2 ( ) ; 

} 


/*********************************************************/ 
/*  SUB-SUBROUTINE  INPUT2()  for  draw2()  */ 

/*********************************************************// 
input2 ( ) 
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{ 

alphal-alphal*180 . O/PI ; 
alpha2-alpha2*180 . O/PI ; 

varl4 ( ) ; 

setcolor (&white) ; 

bxl-460.0;byl-195.0;bx2-638.0?by2«305.0; 
gcx-577 . 0 ygcy-299 . 0 ; 
boxf  ( )  y 

varl3 () ; 

settextclr (iligreen, &blue) ; 
tcx-532 . 0 ; tcy-291. 0  ? 
movtcurabs ( &tcx , & tcy ) ; 
text (" INPUT" ) ; 

settextclr (filiyellow, iblue) ; 
tcx-495 . 0  ? tcy-280 . 0 ; 
movtcurabs ( &tcx, &tcy ) ; 
text ("Encounter") ; 
settextclr (&white, &blue) ; 
tcx«470.0;tcy«270.0; 
movtcurabs ( &tcx , &tcy ) ; 
text ("l.R_det  «"); 
settextclr (filiyellow, &blue) ; 
tcx-495. 0;tcy»2 60.0; 
movtcurabs ( &tcx , &tcy ) ; 
text ( "Warhead" ) ; 
settextclr (Swhite, &blue) ; 
tcx-470 . 0 ; tcy-250 . 0 ; 
movtcurabs (&tcx,& tcy) ; 
text("2.alphal  -")  ; 
tcy-240.0; 

movtcurabs ( &tcx , &tcy ) ; 
text ("3 .alpha2  -") ; 
tcy-2-O.Oy 

movtcurabs ( &tcx, &tcy) ; 
text ( " 4 . N ( frag) «" ) ; 
tcy-220.0; 

movtcurabs ( &tcx , &tcy) ; 
text ("5.  Vo  «"); 
settextclr (fcliyellow, &blue) ; 
tcx-495. 0;tcy-210. 0; 
movtcurabs ( &tcx , &tcy ) ; 
text ("Target") ; 
settextclr ( & white , fcblue) ; 
tcx-470. 0;tcy-200.0; 
movtcurabs ( &tcx , itcy ) ; 
text ("6. length  -") ; 

) 
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/**************************************** *****************/ 
/*  SUB-SUBROUTINE  INPUT21()  for  draw2()  */ 

/*********************************************************/ 

input21() 

{ 

sprintf (s, M%-5.1fH,r_det) ; 
tcx“555 . 0 ; tcy-270 . 0 ; 
varlO ( ) ; 

sprintf (s, w%-6. If M , alphal) ; 
tcy*250. 0; 
varlO ( ) ; 

sprintf (s, " %-6 .  If  " , alpha2 ) ; 

tcy-240.0; 

varlO ( ) ? 

sprintf (s,"%-7. If' ,num) ; 
tcy-230. 0; 
varlO  ()  ; 

sprintf (s,"%-7. If", vo) ; 

tcy-220.0; 

varlO () ; 

sprintf (s, "%-5. If" , length) ; 

tcy-200.0; 

varlO () ; 

} 


/********************************************************/ 

/*  SUB-SUBROUTINE  L0C_2 ( )  for  draw2 ( )  */ 

/********************************************************/ 

loC_2 ( ) 

{ 

var7 ( ) ; 
msl ( ) ; 

} 


/********■*************************************************/ 
/*  SUB-SUBROUTINE  VECTOR2()  for  draw2 ( )  */ 

/*********************************************************/ 


vector2 ( ) 

{ 

alphal-alphal*PI/180 . 0 ; 
alpha2-alpha2  *PI/180 . 0 ; 
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gcx*imx+80 . 0 ;gcy»imy+60 . 0 ; 
movabs ( &gcx, &gcy) ; 

/********  vector  for  Vo  due  to  alphal  *********/ 
setcolor (Sligreen) ; 

dx*dist_l* (vo/vm) *cos (theata+alphal) ; 
dy-dist_l* (vo/vm) *sin (theata+alphal) ; 
dx*dx*l. 439 ; 
lnrel ( &dx, &dy) ; 

/***********  vector  for  Vm  ***********/ 

setcolor(&licyan) ; 

dx-dist_l*cos(theata)  ; 

dy*dist_l*sin(theata) ; 

var9 ( ) ; 

/************  vector  for  Vi  ***********/ 

setcolor (Slipurple) ; 

dx-dist_l*l. 439* (cos (theata)+ (vo/vm) *cos (theata+alphal) ) ; 

dy«dist_l* (sin (theata)+ (vo/vm) *sin (theata+alphal) ) ; 
lnrel (&dx,&dy) ; 

/************  vector  for  vt  ***********/ 

setcolor (fcwhite) ; 

dx«-dist_l*vt/vm; 

dy»0 . 0 ; 

dx=dx* 1.439 ; 

lnrel (&dx, &dy) ; 

/*************  vector  for  Vft  ************/ 

setcolor (Slired) ; 

inqgcur(&gcx, &gcy, &lired) ; 

gcxl*gcx ; gey l=gcy ; 

gcx-imx+80 . 0 ;gcy*imy+60 . 0  ? 

lnabs(&gcx,&gcy)  ; 

> 


/*********************************************************/ 
/*  SUB-SUBROUTINE  VECT0R210  for  draw2 ( )  */ 

/*********************************************************/ 

vector2 1 ( ) 

/**********  vector  for  Vo  due  to  alpha2  **********/ 
setcolor  (Stligreen)  ; 

dx«dist_l* (vo/vm) *cos (theata+alpha2 ) ; 
dy»dist_l* (vo/vm) *sin(theata+alpha2) ; 


dx~dx*1.439; 
lnrel ( &dx , &dy ) ; 

y*************  vector  for*  Vm  ****************y 

setcolor (Slicyan) ; 

dx-dist_l*cos (theata) ; 

dy»dist_l* sin (theata) ; 

var9 ( ) ; 

y ***************  vector  for  Vi  *****************y 

setcolor (Slipurple) ; 

dx»dist_l*l. 439* (cos (theata) +(vo/vm) *cos (theata+alpha2) ) ; 
dy»dist_l* (sin (theata) +(vo/vm) *sin(theata+alpha2) ) ; 
lnrel (&dx, &dy) ; 

y ****************  vector  for  Vt  ***************y 

setcolor (&white) ; 

dx*-dist_l*vt/vm ; 

dy**0 . 0 ; 

dx=dx*l . 439 ; 

lnrel (&dx, &dy) ; 

y ****************  vector  for  vft  ***************y 

setcolor (Slired) ; 

inqgcur (&gcx, &gcy, &lired)  ; 

gcx2=gcx ; gcy2=gcy ; 

gcx=imx+80 . 0 ;gcy=imy+60 . 0 ; 

lnabs(&gcx, &gcy) ; 

y****  triangle  for  fragment  spray  zone  *****y 

movabs ( &gcx2 , &gcy2 ) ; 

lnabs (Sgcxl, Sgcyl) ; 

settextclr (&lired, fiblack) ; 

tcx=«15 . 0  ;  tcy=*3  7 . 0 ; 

movtcurabs (&tcx, &tcy) ; 

text ("Triangle  with  red  lines  shows  fragment  spray  zone") ; 
tcy«25.0; 

movtcurabs ( &tcx , & tcy ) ; 

text ("with  respect  to  target 1"); 

} 


y*********************************************************y 
/*  SUB-SUBROUTINE  INDEX2()  for  draw2 ( )  V 

y*********************************************************y 

index2 ( ) 

{ 

dx=22 . 0 ;dy=0 . 0 ; 
gcx»2  0.0; gey- 3  00.0; 
tcx-47 . 0 ; tcy»2  9  6.0; 
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gcy-290.0; 
tcy*286.0; 
in_vm()  ? 

gcy*280 . 0  ?tcy*276 . 0 ; 
in_vi ( )  ? 

gcy-270.0;tcy-266.0; 
in_vt ( ) ; 

gcy-260 . 0 ; tcy=256 . 0 ; 
in_vft ( ) ; 

} 


y*********************************************************y 
/*  SUB-SUBROUTINE  OUTPUT2Q  for  draw2()  */ 

y*********************************************************y 

output2 ( ) 

{ 

setcolor (&white) ; 

bxl=460 . 0 ; by 1*47 . 0  ?bx2=638 . 0 ;by2=173 . 0 ; 
gcx*577 . o ;gcy*l72 . 0 ; 
boxf  ( )  ; 

varl3  ()  ; 

settextclr (&ligreen, Sblue) ; 
tcx*53  0 . 0 ; tcy*160 . 0 ; 
movtcurabs ( &tcx, &tcy) ; 
text ("OUTPUT") ; 
settextclr (Swhite, &blue) ; 
tcx*470 . 0 ; tcy*150 . 0 ; 
movtcurabs ( &tcx, &tcy) ; 
text ("l.Vftxl  *"); 
tcy*140.0; 

movtcurabs ( &tcx, &tcy) ; 
text ("2 .Vftyl  -"); 
tcy-130. 0; 

movtcurabs ( &tcx , &tcy ) ; 
text("3.Vftx2  *"); 
tcy-120.0; 

movtcurabs (&tcx, &tcy) ; 
text("4.Vfty2  ; 

tcy-110.0; 

movtcurabs (&tcx, &tcy) ; 
text ("S.Vftl  -"); 
tcy»100.0; 

movtcurabs ( &tcx, &tcy) ; 


text ("6.Vft2 
tcy-90.0; 

movtcurabs ( &tcx, &tcy) ; 
text  ('*  7.  Gamma  1  ■")  ; 
tcy-80.0; 

movtcurabs ( &tcx, &tcy) ; 
text ("8. Gamma 2  -") ; 
tcy-70. 0; 

movtcurabs ( &tcx, &tcy) ; 
text("9.Phione  ■") ; 
tcy«60.0; 

movtcurabs ( &tcx, &tcy) ; 
text ( " 10 . Phi two*" ) ; 
tcy-50.0; 

movtcurabs (&tcx, &tcy) ? 
text ("ll.Rho  =") ; 

> 


/*********************************************************/ 
/*  SUB-SUBROUTINE  0UTPUT21()  for  draw2()  */ 

/*********************************************************/ 


output2 1 ( ) 

{ 

gammal=gammal*180. 0/PI; 
gamma2~gamma2*180 . 0/PI ; 
phione-phione* 18  0 . 0/PI ; 
phitwo-phitwo*180 . 0/PI ; 

sprintf (s, "%-8.2f" ,vftxl) ; 
tCX^556 . 0 ;tcy=150 . 0  ? 
varlO ( ) ; 

sprintf (s, "%-8.2f",vftyl) ; 

tcy-140.0; 

varlO ( ) ; 

sprintf (s, "%-8 . 2f" ,vftx2) ; 

tcy«130.0; 

varlO () ; 

sprintf (s," %-8. 2 f" ,vfty2) ; 

tcy-120.0; 

varlO ( ) ; 

sprintf ( s , "%-8 . 2  f " , vf 1 1 ) ; 

tcy«110.0; 

varlO () ; 

sprintf (s,»%-8.2f",vft2)  ; 
tcy-100. 0; 
varlO ( ) ; 

sprintf (s,"%-6. If", gammal) ; 
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tcy-90.0; 
varlO ( ) ; 

sprintf (s,"%-6.1f" ,gamma2) ; 
tcy-80. 0; 
varlO ( ) ; 

sprintf (s , "%-6. If", phione) ; 
tcy»70. 0; 
varlO ( ) ; 

sprintf (s,”%-6. If” ,phitwo) ; 
tcy-60. 0; 
varlO () ; 

sprintf (s,”%-7.4f”,rho) ; 
tcy=50. 0; 
varlO  ()  ; 

gwrite(file2) ; 
bottom () ; 
varl2 ( ) ; 

} 


/*  SUB-SUBROUTINE  LOC_3()  for  draw3()  */ 

loc_3 ( ) 

{ 

trans_factor=0 .02*1.2; 
var8 ( ) ; 
acftl ( )  ; 

} 


/*********************************************************/ 
/*  SUB-SUBROUTINE  VECTOR3()  for  draw3 ( )  */ 

/*********************************************************/ 

vector3 ( ) 

{ 

/****  vector  for  Vft  from  tail  section  of  aircraft  ****/ 

dx=0 . 0 ; dy=0 . 0 ; 

movrel (&dx, &dy) ; 

dx«-trans_factor*vftxl*l . 439 ; 

dyat-trans_f  actor  *vftyl ; 

lnrel (&dx, &dy) ; 

dx*-dx ; dy=-dy ; 

movrel (&dx, &dy) ; 

dx— trans  factor*vftx2* 1.439; 
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dy— trans_f  actor*vf  ty2 ; 
lnrel (&dx, Sidy) ; 
dx— dx ;  dy=-dy ; 
movrel ( Sidx, Sidy) ; 

/****  vector  for  Vft  from  pitot  tube  of  aircraft  ****/ 
dx-8  4 . 0 ; dy=0 . 0 ; 
if  (ht==0 . 0) 
dx=dx ; 
else 

dx=>-dx; 

movrel ( &dx, &dy) ; 

dx=«-trans_factor*vftxl*l.  439  ? 

dy*-trans_f actor*vf tyl ; 

lnrel ( &dx, &dy) ; 

dx=-dx ; dy=-dy ; 

movrel ( &dx, &dy) ; 

dx*-trans_f actor *vftx2*l. 4  39  ? 

dy=-trans_f actor*vfty2 ; 

lnrel (&dx, &dy) ? 

mdist()  ; 

vmt_vector()  ; 

msl_pt()  ; 

} 


/*  SUB-SUBROUTINE  INDEX3 ( )  for  draw3 ( )  */ 

/******************************************★**********■****/ 

index3 ( ) 

{ 

varl3 () ; 
dx*22.0;dy=0.0; 
gcx=52  0 . 0 ; gcy=7  0.0; 
tcx=547 . 0 ; tcy=66 . 0 ; 
in_vft ( ) ; 

gcy=6  0 . 0 ; tcy=5  6 . 0 ; 
in_miss ( ) ; 

gcy=50 . 0 ; tcy=4  6.0; 
in_vmt ( ) ; 

} 


/*  SUB-SUBROUTINE  OUTPUTS ()  for  draw3 ( )  */ 


output3 ( ) 

{ 

setcolor (iwhite) ; 

bxl-520 . 0 ;byl-110 . 0 ;bx2-63  9 . 0 ;by2-2  00 . 0 ; 
gcx-577 . 0 ;gcy-l99 . 0  ? 
box£  ( ) ; 

settextclr(&ligreen,&blue) ; 
tcx-540 . 0  ; tcy-185 . 0 ; 
movtcurabs ( &tcx , &tcy ) ; 
text("Tgt  length"); 

tcy-165.0; 

movtcurabs (&tcx, &tcy) ; 
text ("Hiss  Diet") ; 

tcy-145.0; 

movtcurabs ( &tcx, &tcy) ; 
text ("Gammal") ; 

tcy-125.0; 

movtcurabs  ( &tcx,  &tcy) ; 
text ( ”Gamma2 " ) ; 

> 


/ft********************************************************/ 

/*  SUB-SUBROUTINE  0UTPUT31()  for  draw3()  */ 

/ft********************************************************/ 

output31() 

{ 

settextclr (Swhite, &blue) ; 
sprintf ( s , "%-5 . If" , length) ; 
tcx-560 . 0 ; tcy=17  5.0; 
varl0(); 

sprintf (s,"%-7.3f",m_dist) ; 
tcy-155 . 0 ; 
varlO ( ) ; 

sprintf (s,"%-6. if" ,gammal) ; 

tcy-135.0; 

varlO ( ) ; 

sprintf (s,"% -6. If", gamma2) ; 

tcy-115.0; 

varlO  () ; 

) 
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/a#*********##********************************************/ 
/*  SUB-SUBROUTINE  COORD_3 ( )  for  draw3()  */ 

/*********************************************************/ 

coord_3 ( ) 

{ 

gcx-20 . 0  ;gcy-40 . o  ? 
movabs ( &gcx, &gcy) ; 
gcx-440.0; 
lnabs  ( &gcx ,  &gcy )  ; 
gcx-20 . 0 ;gcy«40 . 0 ; 
movabs ( &gcx, &gcy) ; 
gcy«272.0; 
lnabs (&gcx, &gcy) ; 

for  ( count-0 ; count< 5 ; count++ ) 

{ 

gcx-gcx+8 4.0; gcy-4  2.5; 
movabs (&gcx,&gcy) ; 
gcy-40.01; 
lnaJas  (&gcx,  &gcy)  ; 

} 

gcy-40.0; 

for  ( count-0 ; count<4 ; count++ ) 

{ 

gcx»2  0.01; gcy-gcy+8  4/1.439; 
movabs (&gcx,&gcy) ; 
gcx-24 . 0 ; 
lnabs ( &gcx , &gcy ) ; 

} 

settextclr (Sligreen) ; 
tcx-450 . 0 ;tcy-38 . 0 ; 
movtcurabs ( &tcx , &tcy) ; 
text ("x  [ft]"); 
tcx-3 . 0 ;tcy-280 . 0 ; 
movtcurabs ( &tcx , &tcy ) ; 
text("y  [ft]"); 
tcx-200 . 0 ; tcy-27 . 0 ; 
movtcurabs ( & tcx , &tcy ) ; 
text ( " [ scale-tgt  length ] " ) ; 
deltcur ( ) ; 

) 


/*********************************************************/ 

/*  SUB-SUBROUTINE  REPEAT ()  for  draw3()  */ 

y'*********************************************************// 


repeat ( ) 


{ 

gwrite ( f ile3 ) ; 
iter_l : 
bottom ( ) ; 
bottoml ( ) ; 
iter_3 : 
c-getch ( ) ? 
if  (c— •£') 

{ 

iter_2 : 

setcolor (Sblack) ; 
clr(); 

gread(file2)  ; 
bottom () ; 
bottoml ( ) ; 
iter_4 : 
c»getch ( ) ? 
if  (c—'f) 

{ 

setcolor (Sblack) ; 
clr  ()  ; 

gread(filel) ; 
bottom ( ) ; 
c«getch() ; 
goto  iter_2? 

} 

else 

{ 

setcolor (Sblack) ; 
clr()  ? 

gread(file3) ; 
goto  iter_l; 

} 

} 

else 

{ 

setcolor (&black) ; 

clr(); 

close ( ) ; 

) 


/*********************************************************// 
/*****  MISCELLANEOUS  SUBROUTINE  FUNCTION  *****/ 

</*********************************************************// 


/********** 


aircraft  function 


**********/ 


acftl() 

{ 

varl4 ( ) ; 

setcolor (&white) ; 
gcx-itx+80 . 0 ;gcy=ity+60 . 0 ; 
movabs ( &gcx , &gcy ) ; 
setcolor (filiyellow) ; 

dx*3.0;dy«0.0; 
var2 ( ) ; 

dx»-6 . 0 ; dy-0 . 0 ; 
var2 ( ) ; 

dx»3.0;dy*0.0; 
var2 ( ) ; 

dx»0 . 0 ; dy=-3 . 0 ; 
var2 ( ) ; 

dx*0.0;dy*6.0; 
var2 ( ) ; 

dx«0.0;dy*-3.0; 
var2 ( ) ? 

setcolor (Swhite) ; 
if  (ht=»0.0) 

{ 

gcx*itx+3  6.0; gcy=ity+56 . 0 ; 
} 

else 

{ 

gcx»itx+124 . 0 ;gcy=ity+56 . 0 
} 

movabs ( &gcx , &gcy ) ; 
dx*0.0;dy*3.0; 
var2 ( ) ; 

dx»2 . 0 ; dy=0 . 0 ; 
var2 ( ) ; 

dx«2.0;dy«16.0; 
var2 ( ) ; 

dx-6.0;dy*0.0; 
var2 ( ) ; 

dx»4 . 0 ;  dy=»-12 . 0 ; 
var2 ( ) ; 

dx=12 . 0 ; dy»-l . 0 ; 
var2 ( ) ; 

dx*14 . 0 ; dy*3 . 0 ; 
var2 ( ) ? 

dx-10.0;dy»4.0; 
var2 ( ) ; 

dx-4.0;dy-0.0; 
var2 ( ) ; 

dx-10 . 0  ;dy— 6 . 0 ; 
var2 ( ) ; 


uwi  nnnmmumnuiwn 


dx*8 . 0  ;dy—  1. 5 ; 
var2()  ; 

dx«8 . 0 ;  dy—  2 . 5 ; 
var2 ( ) ; 

dx-4.0;dy*0.0; 
var2 ( ) ; 

dx*0 . 0 ;dy— 0 . 5 ; 
var2 ( ) ? 

dx— 4 . 0 ;dy=0 . 0 ; 
var2 ( )  ; 

dx— 8 . 0 ;  dy— 2 . 5 ; 
var2 ( ) ? 

dx— 8 . 0 ;  dy*- 1 . 0  ; 
var2 ( ) ; 

dx— 2  4.0;  dy*0 . 0 ; 
var2 ( ) ; 

dx=-4 . o ; dy=-2 . 0 ; 
var2 ( ) ; 

dx— 6.0;dy*0.0; 
var2 ( ) ; 

dx— 10 . 0  ;dy*2 . 0 ; 
var2 ( ) ; 

dx— 18 . 0  ;dy=»l .  0 ; 
var2 ( ) ; 

dx*«-2 . 0  ;dy=0 . 0 ; 
var2 ( ) ; 

setcolor(Slired) ; 
dx— 9 . 0 ;  dy=0 .25; 
var2 ( ) ; 
dx»9 . 0 ; 


var2 ( ) ; 
dx— 7 . 5 ; 
var2 ( ) ; 
dx-7 . 5 ; 
var2 ( ) ; 
dx— 8 . 4 ; 
var2 ( ) ; 
dx*8 . 4 ; 
var2 ( )  ; 
dx— 9. 1; 
var2 ( )  ; 
dx-9 . 1 ; 
var2 ( )  ; 
dx— 9 . 0 ; 
var2 ( )  ; 
dx*9 . 0 ; 
var2 ( ) ; 

} 


/************* 

msl() 


missile  function 


****************/ 


{ 

varl4 ( ) ; 
setcolor (fiwhite) ; 
gcx-imx+80 . 0 ; gcy-imy+60 . 0 ; 
movabs ( &gcx , &gcy ) ; 
msl_point ( ) ; 

gcx-imx+80 . 0 ygcy-imy+60 . 0 ; 
movabs ( &gcx, &gcy) ; 
msl_in_l_point () ; 

} 

/**************  missile  point  function 
msl_point ( ) 

{ 

dx-10 . 0  y  dy=0 . 0 y 
var6 ( ) y 

dx— 10 . 0 ;  dy-2 . 0 ; 
varl  ( )  ? 

dx— 3  0 . 0 ;  dy-0 . 0 ; 
varl  ( )  y 

dx-0 . 0 ;  dy— 2 . 0 ; 
varl  ( )  ; 

dx»40 . 0  ydy*0 . 0 ; 
var6  ( )  ; 

dx— 10 . 0  ;dy— 2 . 0  ; 
varl  ( )  ; 

dx— 30.0;dy=0.0; 
varl  ( )  ; 

dx-0 . 0 ; dy=2 . 0 ; 
varl  ( )  ; 

setcolor (Slired) ; 
dx— 2  0 . 0 ;  dy*0 . 0 ; 
varl  ( )  ; 

dx-20.0ydy»0.5; 
var6 ( ) ; 

dx— 18 . 0 ;  dy-0 . 0 ; 
varl  ( )  ; 

dx-l8.0;dy— 1.0; 
var6  ( )  ; 

dx— 18 . 0 ;  dy»o .  0 ; 
varl ( ) ; 

dx-18.0;dy-1.5; 
var6  ( )  ; 

dx— 16 . 0 ;  dy-0 . 0 ; 
varl()  ; 

dx-16 . 0  ydy— 2 . 0 ; 
var6  ( )  ; 

dx— 16.0;dy«0.0; 
varl  ( )  ; 

dx-16.0;dy-2.5; 


var6 ( ) ; 

dx— 14.0;dy-0.0; 
van  ( )  ; 

dx-14 . 0 ;  dy— 3 . 0 ; 
var6 ( ) ; 

dx— 14 . 0 ; dy-0 . 0 ; 
varl ( ) ; 

} 

/************  missile  initial  point  function  ***********/ 
msl_in_ljpoint ( ) 

{ 

setcolor (&cyan) ; 
dx— 3 . 0 ;  dy-0 . 0 ; 
lnrel (&dx,&dy) ; 
dx«6 . 0 ; dy-0 . 0 ; 
lnrel (&dx, &dy) ; 
dx— 3 . 0 ;dy=0 . 0 ; 
lnrel (&dx, &dy) ; 
dx-0 . 0 ;dy=-3 . 0 ; 
lnrel (&dx, &dy) ; 
dx-0 . 0 ;  dy-6 . 0 ; 
lnrel (&dx, &dy) ; 
dx»0. 0;dy— 3.0; 
lnrel (&dx, &dy) ; 
dx»50.0;dy*0.0; 
lnrel (&dx,&dy) ; 
dx— 50 . 0 ; dy=0 . 0 ; 
lnrel (&dx, &dy) ; 

inqgcur (&gcx, &gcy, &cyan) ; 
gcxl-gcx ; gey l=gcy ; 

dx-50.0*cos(theata) *1.439; 
dy-50. 0*sin(theata) ; 
lnrel (&dx, Sidy) ; 
movabs(&gcxl, Sgcyl) ; 

radius-20.0; 
angl-0 . 0 ; ang2-theata ; 
arc ( Sradius , Sangl , &ang2 ) ; 
movabs ( Sgcxl , &gcy 1 ) ; 


/***********  miss  distance  vector  function  *************/ 
mdist ( ) 

{ 

setcolor (Sliblue) ; 
gcx-itx+80 . 0 ;gcy«ity+60 . 0 ; 


movabs ( &gcx , &gcy ) ; 
dx— 84 . 0*minsx/length; 
dy— (84 . 0*minsy/ length) /I . 439 ; 
lnrel (&dx,&dy) ; 

) 


/*********  missile  velocity  vector  function  **********/ 
vmt_vector ( ) 

{ 

setcolor (filiyellow) ; 
dx=vmtx*2*trans_factor*1.439 ; 
dy-vmty*2*trans_f actor; 
lnrel (&dx, &dy) ; 

dx— 2 . 0*dx; 
dy— 2 . 0*dy; 
lnrel (&dx, &dy) ; 

} 


/*********  variable  missile  point  function  *********/ 
msljpt ( ) 

{ 

varl4 () ; 

setcolor ( &white) ; 
inqgcur ( &gcx , &gcy , Swhit e ) ; 
gcxl=gcx ; gey l=gcy ; 
msl_point ( ) ; 
movabs ( &gcxl , &gcy 1 ) ? 
msl_in_ljpoint ( ) ; 

} 


/*************  index  of  Vft  **************/ 

in_vf t ( ) 

{ 

setcolor (&lired) ; 
movabs ( &gcx , &gcy ) ; 
lnrel (&dx, &dy) ; 
settextclr (Slired, Sblack) ; 
movtcurabs(&tcx, &tcy) ; 
text ( HVf t" ) ; 

} 


/*************  index  of  vmt  *************/ 

in_vmt ( ) 

{ 

setcolor (filiyel low) ; 
movabs ( &gcx, &gcy ) ; 
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lnrel (&dx,fcdy) ; 
settextclr (&liyellow, Sblack) ; 
movtcurabs ( &tcx, &tcy ) ; 
text ("Vmt") ? 

) 


/**************  index  of  miss  distance  ***************/ 
in_miss ( ) 

{ 

setcolor(&liblue) ; 
movabs (&gcx, &gcy) ; 
lnrel ( &dx, &dy ) ; 
settextclr ( Sliblue , Sblack) ; 
movtcurabs ( &tcx, &tcy) ; 
text ("Miss  Dist") ; 

} 


/**************  index  of  vi  **************/ 

in_vi ( ) 

{ 

setcolor (Slipurple) ; 
movabs (&gcx, &gcy) ; 
lnrel (&dx,&dy) ; 
settextclr (&lipurple, &black) ; 
movtcurabs ( &tcx , &tcy ) ; 
text ("Vi") ; 

) 


/************  index  of  missile  velocity  ************/ 
in_vm() 

( 

setcolor (Slicyan) ; 
movabs ( &gcx, &gcy) ; 
lnrel (&dx,&dy) ; 
settextclr (Slicyan, Sblack) ; 
movtcurabs ( &tcx,  &tcy ) ; 
text ( "Vm" ) ; 

} 


/*************  index  of  target  velocity  ************/ 
in_vt()  7 

{ 

setcolor (fiwhite) ; 
movabs ( &gcx, &gcy ) ? 
lnrel (&dx, &dy) ; 
settextclr (fcwhite, Sblack) ; 
movtcurabs ( &tcx , &tcy ) ; 
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text ("Vt") ; 
} 


/*************  index  of  fragment  velocity  ***********/ 
in_vo ( ) 

{ 

setcolor (Sligreen) ; 
movabs ( &gcx, &gcy) ; 
lnrel ( &dx, &dy ) ; 
settextclr (fcligreen, Sblack) ; 
movtcurabs ( &tcx , &tcy ) ; 
text (MVoM) ; 

} 

/***************  top  marking  **************/ 

top() 

{ 

h»2 ;w*2 ; 

settext ( &h , &w , &p , &m) ; 
settextclr (Slicyan, Sblack) ; 
tcx*20.0? 
tcy-320.0; 

movtcurabs ( &tcx , &tcy ) ; 

} 


/***************  bottom  marking  ****************/ 

bottom ( ) 

{ 

width  -  13; 
setlnwidth(&width) ; 
setcolor (Sliblue) ; 
gcx*10 . 0 , gcy»10 . 0 ; 
movabs ( &gcx, &gcy ) ; 
gcx-635.0; 
lnabs ( &gcx, &gcy) ; 
varl3 ( ) ; 

settextclr  (Swhite,  Siliblue) ; 
tcx»4 5 0 . 0 ; tcy»7 . 0 ; 
movtcurabs ( &tcx, &tcy) ; 
deltcur() ; 

text ("Hit  any  key  to  proceed") ; 

} 


bottoml ( ) 

{ 

varl3 ( ) ; 

settextclr(&white,&liblue) ; 
tcx«150 . 0 ; tcy«7 . 0 ; 


movtcurabs  ( &tcx ,  &tcy) 
text ("Hit") ; 
tcx-198 . 0 ; 

movtcurabs ( &tcx, &tcy) ; 
text ("to  see  previous  display."); 
settextclr(&lired,&liblue) ; 
tcx-182 . 0 ; 

movtcurabs (&tcx, &tcy) ; 
deltcur ( ) ; 
text("F") ; 


/**************** 
boxf () 


routine  box  function  **************/ 


box ( &bxl , Sibyl , &bx2 , &by2 ) ; 
movabs ( &gcx, &gcy) ; 
flood (&blue) ; 

} 


/**************** 
close  () 


closing  screen 


*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * / 


h=*7  ;w=5  ;p=0  ;m=0 ; 
settext ( &h , &w , &p , &m) ; 
settextclr (Slicyan, Sblack) ; 
tcx=240. 0;tcy=140. 0; 
movtcurabs ( &tcx, &tcy ) ; 
deltcur ( ) ; 
text ("END") ; 
varl2 ( ) ; 

) 


/*********************************************************/ 
/***  routine  varing  functions  ***/ 

/*********************************************************^ 

varl  ( ) 

{ 

if ( ( (3 . 14<theata) && (theata<4 .71) ) | | ( (1.57<theata) && (theata< 
3.14))) 

{ 

dx-dx*cos (theata) -dy*s in (theata) ; 
dy«dx*sin (theata) +dy*cos (theata) ; 
dy— dy; 

lnrel(&dx,&dy) ; 


el 

if ( ( (0.0<theata) &&(theata<1.57) ) | | ( (4.71<theata) &&(theata<6 
.28))) 

{ 

dx»dx*cos (theata) -dy*sin(theata) ; 
dy»dx*sin (theata) +dy*cos (theata) ; 
lnrel ( &dx, &dy ) ; 

} 

) 

var2 ( ) 

{ 

if  (ht==0. 0) 
dx=dx; 
else 
dx— dx ; 

lnrel ( &dx, &dy) ; 

} 

var3 ( ) 

{ 

trans_f actor  =0.2  ; 
if ( (tx<mx) && (ty<my) ) 

{ 

itx=0 . 0 ; ity=0 . 0 ; var5 ( ) ; 

} 

else  if ( (tx<mx) && (ty>my) ) 

{ 

itx=o . 0 ; ity=200 . 0 ; var5 ( ) ; 

} 

else  if ( (tx>mx) && (ty>my) ) 

{ 

imx=0 . 0 ; imy=0 . 0 ; var4 ( ) ; 

} 

else  if ( (tx>mx) && (ty<my) ) 

{ 

imx=0 . 0 ; imy=200 . 0 ; var4 ( ) ; 

} 

) 

var4 ( ) 

{ 

itx=imx+ (tx-mx) * trans_f actor *1 . 439 ; 
ity»imy+(ty-my) *trans_factor ; 

} 

var5 ( ) 

{ 

imx«itx+ (mx-tx) *trans_factor*l. 439 ; 
imy=ity+ (my-ty ) *trans_f actor ; 
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} 


var6 ( ) 

{ 

if ( ( (3.14<theata)&&(theata<4.71) ) | | ( (1. 57<theata) && (theata< 
3.14))) 

{ 

dx«dx*cos (theata) -dy*sin(theata) ; 
dy»dx*s in (theata) +dy*cos (theata) ? 
dy»-dy ; 

movrel ( &dx, &dy) ; 

} 

else 

if ( ( (0. 0<theata) && (theataci. 57) ) | | ( (4 . 71<theata) &&(theata<6 
•28))) 

{ 

dx=dx*cos (theata) -dy* sin (theata) ; 
dy~dx* sin (theata) +dy*cos (theata) ; 
movrel ( &dx, &dy) ; 

> 


) 


var7 ( ) 

{ 

if  ( ( (tx<mx)  &&  (ty<my) )  && (alphaKO .  0) ) 


{ 


imx=280 . 0 ; imy=40 . 0 ; 

} 

else  if  ( ( (tx<mx) &&(ty<my) ) && (alphal>0 . 0) ) 
( 

imx»12  0 . 0 ; imy=18  0 . 0 ; 

} 

else  if  ( ( (tx<mx) &&(ty>my) ) && (alphaKO.  0) ) 
{ 

imx=*o . 0 ; imy=10 . 0 ; 

} 

else  if  ( ( (tx<mx) && (ty>my) ) && (alphal>0 .0) ) 
{ 

imx«280 . 0 ;  imy=*110 . 0  ; 

) 

else  if  (( (tx>mx)&&(ty>my)  )&& (alphaKO. 0) ) 
( 

imx»l0 . 0 ; imy«160 . 0 ; 

) 

else  if  ( ( (tx>mx) && (ty>my) ) && (alphal>0, 0) ) 
{ 

imx«140 . 0 ; imy=10 . 0 ; 

) 

else  if  (( (tx>mx)  &&(ty<my)  )&& (alphaKO.  0) ) 
{ 
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imx«200 . 0 ;  imy**150 . 0 ; 

) 

else  if  ( ( (tx>mx) && (ty<my) ) && (alphal>0 . 0) ) 
{ 

imx«40 . 0 ; imy=40 . 0 ; 

} 


var8 ( ) 

{ 

if  ( ( (txcmx) && (tycmy) ) && (alphal<0. 0) ) 

{ 

itx=12 0 . 0 ; ity=13 0 . 0 ? 

} 

else  if  ( ( (txcmx) && (tycmy) ) && (alphal>0 . 0) ) 
{ 

itx=130 . 0 ; ity=20 . 0 ; 

} 

else  if  ( ( (txcmx) && (ty>my) ) && (alphaKO .  0) ) 
{ 

itx=190 . 0 ;  ity=*170. 0  ; 

} 

else  if  ( ( (tx<mx) && (ty>my) ) && (alphal>0. 0) ) 
t 

itx=120 . 0 ; ity=8  0.0; 

} 

else  if  (( (tx>mx) && (ty>my)  )&& (alphaKO .  0) ) 
{ 

i tx=2  60.0; i ty=7  0.0; 

} 

else  if  ( ( (tx>mx) && (ty>roy) ) && (alphal>0. 0) ) 
{ 

itx=220 . 0 ; ity*150 . 0 ; 

} 

else  if  ( ( (tx>mx) && (ty<my) ) && (alphalcO. 0) ) 
{ 

itx»110 . 0 ; ity*30 . 0 ; 

} 

else  if  ( ( (tx>mx) && (ty<roy) ) && (alphal>0. 0) ) 
{ 

itx*110 . 0 ;  ity=*160. 0  ; 

} 


var9 ( ) 

{ 

dx-dx*1.439  ;• 
lnrel (&dx, &dy) ; 
gcx«imx+80 . 0 ;gcy=imy+60 . 0 ; 
movabs ( &gcx, &gcy) ; 

} 
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varlO() 

{ 

movtcurabs (&tcx, &tcy) ; 
text(s) ; 
deltcur() ; 

} 

varl2 ( ) 

{ 

cratch  ( )  ? 
setcolor (fiblack) ; 
clr(); 

} 


varl3 ( ) 

{ 

h-l?w-l;p«0?m-0? 
settext ( &h, &w, &p , &m) ; 
} 

varl4 ( ) 

{ 

width-1 ; 

setlnwidth(&width) ; 

> 


varl5 ( ) 

{ 

if  (  atan(sy/sx)  <  atan(vmty/vmtx)  ) 
{ 

alphal  -  -alphal  ? 
alpha2  -  -alpha2  ? 

} 

else 

{ 

alphal  -  alphal  ; 
alpha2  -  alpha2  ; 

) 

} 


varl6() 

{ 

if  ( ( (vftxKO.O)  && (vftyl>0. 0) )  |  |  ( (vftxl<0. 0) && (vftyl<0. 0) ) ) 
gamma 1-gamma 1+3 . 14 ; 
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gammal-gammal+6 . 28 ; 
if  (gaamal>6.28) 

ganmal-gasmal-6 . 28 ; 

else 

gamma  1 -ganuna  1 ; 

if  ( ( (vftx2<0 . 0) && (vfty2>0. 0) ) | | ( (vftx2<0. 0) && (vfty2<0. 0) ) ) 
gamma2“gamma2+3 . 14 ; 
else 

gamma2«gamma2+6 . 2 8 ; 
if  (gamina2>6.28) 

gamma2-gamma2-6.28; 

else 

gamma  2 “gamma 2 ; 


APPENDIX  C 


USER'S  MANUAL 


The  following  is  the  user's  manual  of  ENDGAME. 


USER'S  MANUAL 


AE  3251/AE  3705 
AIR  DEFENSE  LETHALITY 


A  STUDY 
Of 

the  ENDGAME 
in 

AIR  DEFENSE 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


I 


INTRODUCTION 


This  homework  problem  is  given  to  the  student  to 
present  an  opportunity  for  him/her  to  see  how  the 
vulnerability  of  an  aircraft  can  be  assessed  in  a  missile- 
based  air  defense  situation.  The  computer  simulation 
program  named  ENDGAME  will  be  used  to  assess  the  missile 
miss  distance,  the  fragment  spray  zone  and  the  best  point 
for  detonation. 


II.  PROBLEM  DEFINITION 

You  are  going  to  conduct  an  endgame  assessment  of  a 
generic  attack  aircraft  on  an  attack  mission  to  destroy  a 
ship  or  ground  vehicle.  Each  student  will  be  given  an 
ordnance  package,  target  and  missile  velocity  and  flight 
path.  This  problem  is  purely  for  instructional  purposes 
and  is  not  based  on  any  actual  or  planned  combat  situation. 
The  aircraft  vulnerability,  missile  characteristics, 
ordnance  package,  flight  path  parameter  limits  have  been 
chosen  only  to  provide  guidelines  for  the  class  problem. 


1.  The  initial  positions  of  the  target  and  the  missile 
are  restricted  to  within  100,000  ft  and  must  be  above 
ground.  That  is,  the  initial  x  and  y  positions  of  the 
target  and  the  missile  must  be  realistic.  The  error 
message  which  tells  you  "  UNREAL!  ASSUME  THAT  IT  BE  ABOVE 
GROUND  AND  REAL."  comes  if  you  type  incorrect  or 

unrealistic  input  data.  Just  input  reasonable  data. 


2.  The  relative  difference  in  distance  between  the 
initial  x  and  y  positions  of  the  target  and  the  missile 
must  each  lie  within  1,000  ft.  If  this  requirement  is  not 
met,  no  warning  message  will  appear,  but  the  results  on  the 
screen  may  be  meaningless. 


3.  The  target  or  aircraft  speed  must  be  greater  than  0 
and  less  than  2000  ft/sec  (  or  Mach  2  )  .  Otherwise,  an 
error  message  tells  you  "  UNREAL!  " .  Just  type  again  the 
appropriate  input. 


4.  Assume  straight  and  level  flight  for  the  air 
target.  The  target  heading,  therefore,  must  be  either  0° 
(left  to  right)  or  180°  (right  to  left) . 


72 


■wwiwwiwuiwuimu  wwu 


n 


N 

h 


< 

* 

i 


i 

I 

I 


I 


5.  The  missile  speed  must  be  greater  than  the  target 
speed.  If  not,  the  error  message  "  UNREAL!"  appears.  So, 
use  reasonable  input  data. 


6.  The  input  for  the  elevation  angle  of  the  missile, 
9,  must  also  be  reasonable.  If  the  x  and  y  positions,  for 
example,  of  the  target  is  greater  than  those  of  the 
missile,  and  if  9  is  given  as  145°  or  210°,  then  the  error 
message  will  warn  you  to  use  the  appropriate  elevation 
angle  in  order  to  allow  the  target  to  be  approached. 


III.  EXECUTING  ENDGAME 

You  can  execute  ENDGAME  on  any  IBM  PC/AT  personal 
computer  or  compatible  with  a  monitor  graphic  resolution 
640  x  350.  If  you  have  a  different  system  on  a  personal 
computer,  a  few  modifications  to  the  program  are  required 
in  order  to  run  ENDGAME.  Further  details  describing  these 
modifications  are  given  in  this  user's  manual.  The 
following  instructions  tell  you  how  to  run  ENDGAME.  Three 
different  displays  will  come  out  from  ENDGAME  program. 


1.  Turn  on  the  system  power  of  your  personal  computer. 

2.  Insert  the  diskette  which  contains  the  execution 
file,  the  source  file  and  the  auxiliary  files. 

3.  Run  the  executable  file,  endgame.exe,  by  simply 
typing  "ENDGAME".  Then,  you  will  see  the  opening 
screen. 

4.  Just  follow  the  instructions  shown  on  the  screen  to 
do  the  simulation. 

5.  Do  not  hit  any  key  to  proceed  until  the  message 
line  comes  on  the  bottom  of  the  each  screen.  This 
means  that  the  current  display  is  being  written  to 
the  specified  file  while  the  In-Use  light  on  the 
diskette  driver  comes  on. 

6.  You  will  be  given  the  opportunity  to  see  the 
previous  screen  by  hitting  "F"  key  when  the  best 
point  for  detonation  screen  is  displayed. 


IV.  OUTPUT 


Make  sure  that  your  printer  is  ready  if  you  want  to 
obtain  the  hard  copy  of  the  display.  Just  hit  the  "prt  sc" 
key  while  holding  the  shift  key,  or  hit  the  "print  screen" 
key  by  itself  on  other  special  computer  keyboards. 


V.  WAYS  TO  MODIFY  THE  PROGRAM 

To  modify  the  program  in  order  to  use  the  system 
that  you  have  or  you  are  going  to  use,  you  must  make  the 
following  changes. 

1.  Among  the  external  variables  declared,  the  number 
of  mode  and  page  which  are  declared  by  integer  should  be 
varied  by  the  systems  that  you  have.  That  number  depends 
upon  how  much  memory  is  installed  on  the  board  and  what 
kind  of  graphic  device  driver  you  have.  See  section  7  in 
Function  Description  of  HALO  manual  (Version  2.26). 

2.  The  graphics  device  driver  must  also  be  changed. 
See  the  HALO  manual. 
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